Engineering optical soliton bistability in colloidal media 
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We consider a mixture consisting of two species of spherical nanoparticles dispersed in a liquid 
medium. We show that with an appropriate choice of refractive indices and particle diameters, it 
is possible to observe the phenomenon of optical soliton bistability in two spatial dimensions in a 
broad beam power range. Previously, this possibility was ruled out in the case of a single-species 
colloid. As a particular example, we consider the system of hydrophilic silica particles and gas 
bubbles generated in the process of electrolysis in water. The interaction of two soliton beams can 
lead to switching of the lower branch solitons to the upper branch, and the interaction of solitons 
from different branches is phase-independent and always repulsive. 
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I. INTRODUCTION 

Spatial optical solitons are formed when the change of 
nonlinear refractive index induces an effective lensing ef- 
fect that balances diffraction of the laser beam |l| . When 
a laser beam passes through a colloidal medium com- 
posed of liquid suspension of dielectric nanoparticles, the 
optical gradient force acts against particle diffusion, in- 
creasing the refractive index in the regions of higher light 
intensity. The corresponding local change of the refrac- 
tive index is of the self-focusing type, and allows for cre- 
ation of spatial optical solitons in the form of self-trapped 
beams, as was demonstrated in both theoretical and ex- 
perimental studies 0-0|- We note that optical solitons 
have been also observed in other soft matter systems, 
including liquid crystals |8| and polymers [9j. 

Recently, it was shown theoretically Q that the optical 
response of a colloidal medium in the hard-sphere approx- 
imation allows for the existence of two different stable 
soliton solutions for the same beam power, i.e. the soli- 
ton bistability of the first kind [13: El- Moreover, inter- 
actions of these bistable solitons have interesting proper- 
ties, not found in other optical soliton systems [12 ]. Soli- 
ton switching through collisions and phase-independent 
repulsive interactions were demonstrated. The soliton 
bistability, however, has been only predicted in the one- 
dimensional case [6|, corresponding to a planar waveg- 
uide setup or surface waves. In contrast, it was shown 
that in the two-dimensional case the lower soliton branch 
becomes unstable [13j due to the intrinsic instability of 
multidimensional solitons in Kerr media. Moreover, even 
in the one dimensional case the power range for which 
bistable solitons exist is limited. 

In this paper, we demonstrate that the introduction of 
a second colloidal species allows for the bistability of two- 
dimensional soliton beams. The idea of "engineering" the 
nonlinear response of a colloidal medium by mixing sev- 
eral components was first proposed in [4| . We show that 
with an appropriate choice of refractive indices and parti- 
cle dimensions, it is possible to observe the bistability in 
a broad beam power range. As a particular example, we 



consider the system of hydrophilic silica particles and gas 
bubbles generated in the process of electrolysis in water. 
In analogy with the one-dimensional case, the interaction 
of soliton beams can lead to switching to the upper soli- 
ton branch, and the interaction of solitons from different 
branches is phase-independent and always repulsive. 

II. MODEL 

We consider a mixture consisting of two species of col- 
loidal particles dispersed in a background liquid medium. 
We assume that the first (second) group of particles is 
characterized by a refractive index higher (lower) than 
the background medium index. For the sake of clarity, 
we will refer to the two groups as to "particles" and "bub- 
bles" , although each of the species can have the form of 
solid state, liquid, or gas, as long as the condition im- 
posed on their refractive indices is fulfilled. In particu- 
lar, bubbles can be replaced eg. by a solution of colloidal 
microshells [l4j to enhance the system stability. We also 
assume purely dielectric response of all the components 
to the laser light and negligible absorption. 

The optical gradient force induced by a laser beam will 
have the effect of locally increasing the concentration of 
particles and decreasing the concentration of bubbles, as 
was shown previously If the background concentra- 
tion of bubbles (that is, concentration in the absence of 
laser light) is low enough to describe them as a system 
of noninteracting particles (ideal gas) , the same descrip- 
tion will hold in the presence of an optical beam. On 
the other hand, concentration of particles will increase, 
and interactions between them can become important, 
even if they were negligible at the outset. We take into 
account interactions between the particles in the hard 
sphere approximation. We also assume that the (other- 
wise negligible) interaction between particles and bubbles 
is repulsive, which ensures stability of the system. 

The model of nonlinear laser beam propagation in a 
colloidal suspension in the hard sphere approximation 
was elaborated in Q . In the recent experiment @ it was 
shown that the optical response of colloidal suspensions 
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of polystyrene beads can be substantially different from 
that predicted by this approximation. We notice, how- 
ever, that other types of colloidal systems can be rea- 
sonably well described by the hard-sphere model (l5j . 
Moreover, it is plausible to expect that the phenomenon 
of soliton bistability can be observable also in the case 
of "soft" interaction. The necessary requirement here is 
the saturation of colloidal particle concentration at high 
packing fractions, which must occur due to limited avail- 
able volume. We can therefore treat the hard-sphere ap- 
proximation as the first step to describing and under- 
standing the physics of more complicated systems. 

In our case, the refractive index of colloidal particles 
n p , bubbles rib and the background index riB fulfill the 
condition n p > riB > n b - This assumption is neces- 
sary for observation of the phenomena described below. 
We also assume that the particle and bubble diameters 
are much smaller than the laser wavelength in the back- 
ground medium, dp, <C Xq/ub (Rayleigh regime). The 
osmotic pressure can be calculated from the equation of 
state [iff 

—Z = Z v (r}„), ^ = p,b (1) 

Pv 

where the index v is replaced by p for solid particles 
or b for bubbles, /3 = l/k^T, HV is the osmotic pres- 
sure, p v is the colloidal particle (or bubble) concentra- 
tion, Z v (rj) is the compressibility, and r\ v — p v V v is 
the packing fraction, where V v denotes the particle vol- 
ume. For bubbles we take the ideal gas compressibil- 
ity Zb(r/b) = 1- For solid particles interacting through 
a hard sphere potential, the Carnahan-Starling formula 

^ p (»7p) ~ (i + v P + vl~ r?p)/(i - ^p) 3 s ives a ver y S ood 

approximation up to the fluid-solid transition at rj p w 0.5 
16] . This phenomenological formula is in agreement with 
exact perturbation theory calculations as well as molec- 
ular dynamics simulations. 

We assume for the time being that the gradient of the 
concentration of colloidal particles or bubbles p(r) is lo- 
cally parallel to x, and consider a small box of volume 
dV = dccdS 1 , with length dir. and normal surface dS. The 
difference in the osmotic pressure exerted on the right 
and left surface dll gives rise to an effective force acting 
on the colloidal particles Fi nt . It is equal to the external 
force that is necessary to sustain the concentration gra- 
dient, and dll = -F int /dS = -f int pdV/dS = -f int pdx, 
where f[ nt is the average force acting on a single particle. 
Using Eq. (fTJ) we get d(pZ)/dx — —f- m tp/3- The particle 
current density is equal to 

~3 = PPifcx + 7 int) = PP~7cx - DW(pZ) , (2) 

where p is the particle mobility, and D = p//3 is the 
diffusion constant. In the ideal gas limit, this equation 
becomes Eq. (3) of 

Let m = n u /riB be the ratio of the colloidal particle 
(or bubble) refractive index to the background refractive 
index. Polarizability of a sphere is given by 



where V v — (7r/6)e^ is the sphere volume and 8 V = {m 2 — 

l)/(m 2 +2). If we look for the steady state (j v = 0) in the 

presence of optical field gradient (f C x = (a„/4)V|I?| 2 ) 
we obtain 



,/3d|£| 2 _ d{ Pv Z v ) 



d.r 



da; 



(4) 



which can be solved analytically to give the dependence 
between \E\ 2 and the packing fractions r\ v 



j\E\ J = — [<?(%)- <?(%,)] = 
4 a p a b 



VobJ 



(5) 



where g(j]) = (3— 77) /(l — r/) 3 + lnry, and 770^ is the back- 
ground packing fraction of the v component. 

Assuming relatively low packing fractions, the corre- 
sponding nonlinear refractive index change can be ap- 
proximately calculated using the Maxwell-Garnett for- 
mula [17]. For low refractive index contrast (n„/riB ~ 1) 
we have 



"off = £ eff 



£b ^1 + 3^5„?7 I/ j 



(6) 



Substituting this formula to the Helmholtz equation 
\7 2 E + k 2 n 2 s E — 0, we obtain the propagation equa- 
tion for the slowly varying envelope of electric field u(r) 
defined by E(r) = (2/y/]3)u(r ) exp (ik noz), where n — 



1 



du 



dz 2fcoriQ 



V 2 u + 3fc 2 2J 8 v (t) v - rj 0u )u 



(7) 



where k — fco n B = 27TOb/Ao and the additional last 
term on the left hand side accounts for damping due 
to Rayleigh scattering from the dielectric spheres. The 
damping coefficients are given by j„ = 27r 5 p l/ <5 2 d^/(3A 4 ) 
[l8|, where A = Xq/ub- Additionally, for steady state 
solutions, relation ([5]) gives 



a b \u\ 2 = ln(r7b/r/ob) 



(8) 



at each point of space. 

We renormalize spatial coordinates according to 
(x,y) = (2/3k 2 )^ 2 x (x,y) and z = (2n /3kn B ) x z, 
obtaining 



du 1 
'^ + 2 



^\u + '^8 u {r] v -T] 0v )u+l- ^T^u = , (9) 



where the renormalized damping coefficients are 
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T v = \^ \l + 3$>^ (± 



(10) 



(3) 



From Eq. (JSJ and formula pH|) we conclude that the effect 
of scattering losses depends strongly on the ratio of the 
particle size to the laser wavelength. 



In the following, we consider a particular system 
of hydrophilic silica particles and hydrogen or oxygen 
nanobubbles generated in the process of electrolysis in 
water pjj]. Nanobubbles solutions generated in this way 
can remain stable for several days. The interaction be- 
tween silica particles and bubbles is repulsive indepen- 
dently of the distance [20] , which ensures the system 
stability. We assume that silica particles have diameter 
dp = 50 nm and n p = 1.45 and bubbles have diameter 
c?b = 100 nm and = 1, while the water refractive index 
is iib = 1.33. The background concentrations of particles 
and bubbles in regions of low light intensity are taken as 
^o P = 10 -3 and 770b = 10~ 2 , which is consistent with 
the experiments [19|. In real systems, it is not possible 
to prepare monodisperse solutions. We assume that the 
dispersion of sizes is relatively low, which allows us to 
use the present model with d p and dh equal to average 
or effective sizes. The relatively long laser wavelength 
Ao = 1064 nm allows for lower Rayleigh scattering losses. 

The typical dependence of packing fractions of both 
components on the light intensity as well as the total 
induced refractive index change An = Uefi — no are pre- 
sented in Fig. [TJ The bubble packing fraction is depleted 
almost completely already at relatively low light inten- 
sities due to the large bubble size and the high polariz- 
ability. On the other hand, the concentration of silica 
particles remains low until the intensity reaches a cer- 
tain higher value. Further on, the concentration of the 
particles increases exponentially and finally saturates as 
the packing fraction becomes larger than 10%. In re- 
sult, thanks to the appropriate choice of particle sizes 
and concentrations, the refractive index dependence has 
a form of two "steps" . The first step is a consequence 
of saturation of bubble induced nonlinear index change 
following bubble depletion while the second step is 
an effect of saturation of particle induced nonlinearity at 
high packing fractions due to limited available volume [6j . 
In the next section we show that this artificially prepared 
nonlinearity allows for the bistability of two-dimensional 
soliton beams. 
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FIG. 1: (Color online) (a) Packing fractions of silica particles 
(dashed) and bubbles (dash-dotted), and the total induced re- 
fractive index change (solid) vs. the light intensity. The inset 
presents the dependence in the low intensities range, (b) The 
same presented in a bilogarithmic scale. Thanks to the appro- 
priate choice of particle sizes and concentrations (see text), 
the refractive index dependence has a form of two "steps", 
and supports two-dimensional soliton bistability. See text for 
values of parameters. 



III. TWO-DIMENSIONAL SOLITONS AND 
SOLITON BISTABILITY 

We now consider two-dimensional spatial soliton solu- 
tions of the attenuation- free version of Eq. ©. We look 
for localized solutions in the form of circularly symmetric 
beams u(r) = A{r) exp(iKz) under the condition T v = 0. 
The propagation equation ^ reduces to 



KA+ \ {& + \^) S ^u-Vo»)A = . ( 



11) 



The soliton profiles can be obtained using numerical re- 
laxation methods, 



21|. The dependence of the 
3/r|w| 2 dr/P 

vs. the propagation constant k is displayed in Fig. HJa). 



soliton power P 



see eg 

J |u| 2 dr and width W 



The branches of stable solutions correspond to positive 
slope dP/dn > while unstable solutions are character- 
ized by a negative slope jl0| . The picture shows two sta- 
ble branches, and within the power range P « 4 — 17 W 
stable solutions corresponding to both of the branches 
exist. These bistable solitons fulfill all the three stabil- 
ity conditions required for their robustness during colli- 
sions 22] . 

This picture should be compared with the results ob- 
tained in the case of a single colloidal component. The 
effect of soliton bistability was predicted in the one- 
dimensional case @, corresponding to a planar waveg- 
uide setup or surface waves, however the range of pow- 
ers supporting the bistability was significantly smaller. 
Moreover, it was shown that in two-dimensional case the 
lower soliton branch becomes unstable uM due to the 
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FIG. 2: (Color online) (a) Soliton power (solid) and width 
(dash-dotted) vs. the propagation constant k. Two stable 
branches with dP/dn > exist. The dashed lines depict 
power dependence in the cases when one of the colloidal com- 
ponents (particles or bubbles) is removed from the system. 
Bottom panels show the soliton intensity profiles (solid) and 
colloidal particle (dashed) and bubble (dash-dotted) packing 
fractions for two bistable solitons carrying power P = 10 W 
from (b) the lower stable branch and (c) the upper stable 
branch. In (c), the bubble packing fraction is negligible. No- 
tice the difference in the width scale. 
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FIG. 3: (Color online) Collision of two identical out of phase 
solitons (A</3 = 7r) from the lower branch (left column) and 
of two solitons from different branches (right column). The 
power in each of the beams is equal to P — 15 W (left) and 
P = 12 W (right), and the angle between the beams is 5° 
and 2.5°, respectively. The consecutive frames correspond to 
propagation distances z — (top row) , z = 165 (im (mid- 
dle row), and z = 425 ptm (bottom row), In the first case, 
repulsive interaction leads to switching of the solitons from 
the lower (a) to the upper branch (c). In the second case, 
the interaction is also repulsive, independently of the relative 
phase between the solitons, and leads to destruction of the 
wider but weaker soliton. 



intrinsic instability of multidimensional solitons in Kerr 
media. The introduction of the second colloidal species is 
therefore necessary for the bistability of two-dimensional 
solitons. Indeed, if one of the components is removed 
from the system, only one of the stable branches remain, 
see dashed lines in Fig.^a). It is clear that the bubbles 
determine the properties of the system for low values of k 
and low light intensities, while the particles are the main 
acting component in the regime of high k. The bistabil- 
ity is the result of combination of the effects that these 
two species have on the nonlinear response of the system. 

In Figs. [H[b,c) we present examples of bistable soliton 
profiles corresponding to the soliton power P = 10 W. 
It is clear that the light intensity of the lower branch 
soliton (b) corresponds to the first "step" from Fig. [TJ 
while the upper branch soliton (c) corresponds to the 
second "step" . Hence, one can call the lower branch soli- 
tons "bubble solitons", while the upper branch solitons 
are "particle solitons" , if referring to the main stabilizing 
component. Since the width of the soliton from the lower 
branch is approximately 7x larger than the width of the 



soliton from the upper branch carrying the same beam 
power, these solitons can be easily distinguished in exper- 
iment. We note that in the case of the lower branch soli- 
tons packing fractions of both components are below 1%, 
while the packing fraction of particles is as high as 60% 
in the center of the soliton beam from the upper branch. 
This value indicates that the simple Carnahan-Starling 
model that we have used for description of particle inter- 
actions breaks down, an d ap pearance of ordered dense 
phase can be anticipated [16| . We note however that the 
breakdown of our model should not lead to qualitative 
changes in the beam propagation and soliton properties, 
since the discrepancy in the particle concentration, which 
determines the nonlinear response, is relatively low. On 
the other hand, the packing fraction of bubbles is always 
below 1% (the background packing fraction), which con- 
firms that our assumption of low bubble concentration is 
well justified. 
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IV. SOLITON INTERACTIONS 

We proceed to the investigation of interactions of the 
solitons from the two bistable branches. We consider a 
collision of two soliton beams angled towards each other 
which initially have the form of two stationary soliton 
solutions, Ui(x, y) and u 2 (x, y), separated by a distance 
2x$ large in comparison to their widths. The solitons 
have imprinted opposite linear phases ko, which resemble 
the initial beam tilt, and a constant phase difference At/? 

u(x, y,z = 0) = Ul (x+x , y)e lkoX +u 2 (x-x ,y)e- tkoX+tAv 

The initial profiles are taken as solutions to the un- 
damped equation (jllj) . and the evolution of the beams 
is modelled using the full equation © with the scatter- 
ing losses included. 

In Fig. [3] we present results for interaction of two iden- 
tical out of phase solitons = 7r) from the lower 
branch (a,b,c) and a collision of solitons from different 
branches (d,e,f). We find that, similarly as in the one- 
dimensional case [Hj], interaction of two solitons from 
lower branch can lead to switching to the upper branch, 
Fig. [3Jc). However, we were not able to observe a sim- 
ilar phenomenon in the interaction of two solitons from 
different branches. Instead, the lower branch soliton is 
destructed in most cases, see Fig. [210 • On the other 
hand, we find that the upper branch soliton is much more 
robust and appears in the same form after the collision. 
Nevertheless, we found that the interaction between soli- 
tons from different branches is phase independent and al- 
ways repulsive, in analogy with the one-dimensional case 
[12| . This kind of interaction has been explained within 



the model of effectively incoherent beams due to the large 
difference in propagation constant k. Despite that each 
of the solitons appears as an attractive potential well for 
the other soliton, repulsive interaction occurs if the col- 
lision angle is small enough [l2|, [23| . To our knowledge, 
we present the first example of soliton repulsion from at- 
tractive potential in two dimensions. 



V. CONCLUSIONS 

We have derived the model equations for description of 
a system of two colloidal components dispersed in a liquid 
medium in the presence of coherent laser light. We have 
shown that with an appropriate choice of refractive in- 
dices and particle dimensions, it is possible to observe the 
phenomenon of optical soliton bistability in two spatial 
dimensions in a broad beam power range. Analogously 
as in the one-dimensional case, the interaction of soliton 
beams can lead to switching to the upper soliton branch, 
and the interaction of solitons from different branches is 
phase-independent and always repulsive. The presented 
results can have implications for the ex per iments on op- 
tical solitons in soft matter media 
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